■fi. 


i  AD 

a 

**■»—»?  -t-wnnn  Mir  Wnr  w  nrrm  niwtrawti  > 


rsr's  'rjr~  o  Tf  ?"  tvt  vr  «"*<  i 

I  * 


u  fa  £j  r  'Uik.  i 


"VT-7 i 02 


3Y  RADIATION  AND  CONVECTION  FROM  ANNULAR 
TRAPEZOIDAL  PROFILE  ON  A  CYLINDER 


by 

GARRY  C.  CAROFANO 

AND 

MANUEL  V.  AVILA 


{  DISTRIBUTION 

STATEMi.ii  J,‘  iT  | 

|  Approved  for 

public  itJease;  } 

|  DirRruvjjic. 

;  UnEtnited 

1 

**»»  •“’*  ““w  •»»*■*■  ■»  >•**  »v4i>  mMMSMiw.’hii  iuetb 

JANUARY  1971 


Ipr  |  A-nfif 

m  Ha  S*  Pp  il  P&  &  1 1 


i|  #,  WR  <m 


LUm 


|HK65<S^fltt  ' 

'ti?  t'  *i?  H  ltt)&  I 

WATERVLIET-r.NEW  ' 


AMC8S  Mo.  5523. 11.50300.01 
ns  PROJECT  HO,  W562803A299 

Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

Springfield,  Va.  221M 

BEST  AVAILABLE  COPY 


»v. 

c-i  0 


-us  troy  this  report  when  it  is  no  longer  needed.  Do  not  return  i. 


L--  ■,  I  ‘-'I  Fi 


fhc  findings  in  this  report  a re  not  to  bs  construed  as  an  officio. 

•.‘C.?;:rnt  of  the  Army  position. 


BEST  AVAILABLE  COPY 


DOCUMENT  CONTROL  DATA  -RAD 

tf  clDooJffeofloff  mt  f Iff#*  k«V  •#  aftfliHl  Ml 


I.  ORICINATIN*  ACTIVITY 

Natervliet  Arsenal 
Watervliet*  N.Y.  12189 


■Ml  feo  nliW  •*«  At  tiwtll  rtitW  It  (ItttllM: 


HEAT  TRANSFER  BY  RADIATION  AND  CONVECTION  FROM  ANNULAR  FINS  OF  TRAPEZOIDAL  PROFILE  ON 
A  CYLINDER 


4.  DitcntPTi  v(  noth  (Tfpm  ml  rtptrt  anti  twhtltt  tiMij 

Technical  Report 


AuTHoniti  (rim  MM,  alMt  initial,  teal  mmj 

Garry  C.  Carofano 
Manuel  V.  Avila 


AT  OATS 

January  1971 


M.  CONTRACT  OR  OAANT  NO. 

AMCMS  No.  5523.11.58300.01 

v  aaojtcr  no. 

DA  Project  No.  1W562603A299 


7  A.  TOTAL  NO.  OP  PAKI  If*.  NO.  OF  N(FI 


•A  ONiaiNATON'l  AIPORT  NUMUPlil 


WVT-71U2 


10.  OIOTMWUTION  ITATIUKNT 

This  document  has  been  approved  for  public 
unlimited. 

release  and  sale;  its  distribution  is 

U.  BRONtORlNC  MILITARY  ACTIVITY 

U.  S.  Army  Weapons  Command 

The  equations  describing  steady  heat  transfer  by  radiation  and  convection 
from  annular  fins  of  trapezoidal  profile  on  a  cylinder  are  derived.  Mutual  irradiation 
between  adjacent  fins  is  considered  and  the  appropriate  differential  angle  factors  are 
presented.  A  constant  convection  coefficient  is  assumed  over  all  surfaces  and  the  fins 
are  treated  as  gray  body  radiators.  The  governing  integral  equations  are  approximated 
by  numerical  quadrature  formulas  and  the  resulting  set  of  non-linear  algebraic  equa¬ 
tions  is  programmed  for  solution,  A  Fortran  program  listing  is  included.  Some  results 
are  obtained  for  the  XM144E1  mortar  tube. 


*  «*»«1473 


MPLACII  DO  PORM  I4TR.  I  JAM  M.  WHICH  1C 
OMOLITI  worn  A  AMT  INC. 


Unclassified 

i*c«iltr  CtRRtincation 


Unclassified 

wity  CUwllleitlw 


Heat  Transfer 
Steady  State 
Fins 

Cylinder 

Conduction 

Convection 

Radiation 


HEAT  TRANSFER  BY  RADIATION  AND  CONVECTION  FROM  ANNULAR 
FINS  OF  TRAPEZOIDAL  PROFILE  ON  A  CYLINDER 


ABSTRACT 

Hie  equations  describing  steady  heat  transfer 
by  radiation  and  convection  from  annular  fins  of 
trapezoidal  profile  on  a  cylinder  are  derived, 
fiitual  irradiation  between  adjacent  fins  is  considered 
and  the  appropriate  differentia},  angle  factors  are 
presented.  A  constant  convection  coefficient  is 
assuned  over  all  surfaces  and  the  fins  are  treated  as 
gray  body  radiators.  The  governing  integral  equations 
are  approximated  by  nunerical  quadrature  formulas  and 
the  resulting  set  of  non-linear  algebraic  equations  is 
prog  rawed  for  solution.  A  Fortran  program  listing  is 
included.  Same  results  are  obtained  for  the  XM144E1 
mortar  tube. 
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A  description  of  the  most  frequently  used  symbols  is  given  below. 
The  remaining  symbols  used  in  the  text  are  defined  when  they  are 
introduced. 


SYMOL 


DESCRIPTION 

Convective  film  coefficient  [Btu/hr-ft2-°R] 

r  2 

Convection  parameter  ■  **  Ri 


■  h  Rj 

km  cos* (6  «•  (Rq  -  Rj)  tan  •<) 


Thermal  conductivity  of  air  [Btu/hr-ft-#R] 

Thermal  conductivity  of  fin  material  [Btu/hr-ft-°R] 
Distance  between  fin  tips  [ft] 


Radiation  parameter 


r  _ >  —•> 

m  Ri  o  Tb 

kn  cos*  (6  ♦  (Ro  -  Ri)  tan^) 

•  m 


Dimensionless  parameter  [  *  I/Rq] 


qr 

*r 

q(r) 

qb(*b) 

Q(^) 

r 

Ri 

s 

So 


Dimensionless  parameter  [  -  Rj/Rq] 

Total  heat  transfer  rate  [Btu/hr] 

Total  heat  transfer  rate  per  unit  length  [Btu/hr-ft] 
Local  fin  radiative  heat  flux  [Btu/hr-ft2] 

Local  base  radiative  heat  flux  [Btu/hr-ft2] 
Dimensionless  variable  [  -  q(r)/ol  ] 

Dimensionless  variable  [  -  q^Cx^/oT^] 

Radial  coordinate  [ft] 

Base  cylinder  radius  [ft] 

Fin  tip  radius  [ft] 

Fin  spacing  [  ■  L  ♦  2  6] 

Fin  spacing  on  XM144E1  mortar  [see  Figure  1] 
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SYMBOL 


DESCRIPTION 


T(r) 

T*(() 

*b 


GREEK  SYMBOLS 
a 
6 

6* 

c 

c 

3 

o 

Tb 

SUBSCRIPTS 

1 

2 

b 

e 


Base  cylinder  temperature  [*R,  unless  otherwise 
noted] 

Environment  temperature  [*R,  unless  otherwise  noted] 
Dimensionless  parameter  [  *  Tg/T^] 

Local  fin  temperature  [  *R  ] 

Dimensionless  variable  [  -  T(r)/Tfc] 

Longitudinal  coordinate  [ft] 

Fin  angle  [degrees] 

Fin  tip  half- thickness  [ft] 

Dimensionless  parameter  [  *  <5 /Rq] 

Emissivity 

Dimensionless  variable  [  -  t/Rq] 

Fin  system  effectiveness 

o  2  4 

Stefan-Boltzmann  constant  [1.714xlO’9Btu/hr-ft  -*R  ] 
Dimensionless  variable  [  ■  15/ Rq] 

Refers  to  fin  1 
Refers  to  fin  2 
Refers  to  base  cylinder 
Refers  to  environment 
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INTRODUCTION 

As  an  tid  in  formulating  the  present  problem  some  of  the  thermal 
operating  characteristics  of  mortars  will  first  briefly  be  reviewed. 

then  a  mortar  is  fired  a  fraction  of  the  chemical  energy  released 
in  the  propellant  combustion  process  is  absorbed  as  thermal  energy  at  the 
bora  surface.  In  a  cold  weapon  a  portion  of  this  energy  is  retained  by 
the  barrel  in  the  form  of  internal  energy;  the  remainder  is  conducted 
through  the  tube  wall  and  dissipated  at  the  outer  surface  by  convective 
aid  radiative  exchange  with  the  surrowdings.  In  a  sustained  firing 
situation  the  energy  storage  results  in  an  increase  in  average  barrel 
temperature,  although,  for  moderate  firing  rates,  a  quasi-steady  radial 
temperature  distribution  eventually  develops  through  a  major  portion  of 
the  tuba  wall;  the  rate  of  energy  input  at  the  bore  surface  is  then  just 
balanced  by  the  dissipation  rate  at  the  outer  surface. 

Also  characteristic  is  the  development  of  an  appreciable  non-uniformity 
of  temperature  in  the  axial  direction  which  is  related  to  the  complicated 
nature  of  the  internal  gas  flow;  the  axial  temperature  distribution  increases 
with  distance  from  the  breech,  attains  a  maximal  and  then  decreases  monot- 
cnically  towards  the  nuzzle.  This  distribution  also  becomes  approximately 
steady  with  time  in  the  outer  regions  of  the  barrel  wall. 

The  design  of  a  mortar  is  controlled  to  a  considerable  extent  by  the 
mechanical  and  thermal  behavior  of  the  tube  segment  located  in  the  vicinity 
of  the  axial  temperature  maximum.  It  is  at  this  point  that  the  tube  strength 
is  at  a  minimuR  because  of  the  inverse  relationship  between  temperature  and 
the  mechanical  properties  of  the  material.  Through  a  suitable  failure 


criterion  the  internal  gas  pressure  and  the  tube  wall-ratio  can  be 
related  to  the  yield  stress  of  the  material.  Since  the  latter  quantity 
has  a  unique  dependence  on  temperature,  the  maximum  sustained  operating 
temperature  which  can  be  permitted  is  that  corresponding  to  the  minimun  t 

yield  stress  which  will  just  maintain  the  structural  integrity  of  the  tube 
under  the  specified  internal  loading  conditions.  The  maximum  heat  dissi¬ 
pation  rate  at  the  outer  surface  is  also  determined  by  this  particular 
value  of  temperature  and,  since  this  is  balanced  by  the  energy  input  at 
the  bore  surface,  a  maximum  firing  rate  is  implied. 

It  would  be  desirable  to  design  a  mortar  with  a  permissable  firing 
rate  which  exceeds  the  maximum  rate  at  which  individuals  could  load  it. 

9 

Although  this  goal  is  far  from  being  attained  by  present  day  mortars, 
significant  advances  have  been  made  by  increasing  the  heat  transfer  rate 
from  the  outer  surface  through  the  additions  of  fins.  An  example  of  such 
a  design  is  the  XM144E1  mortar,  a  section  of  which  is  shown  in  Figure  1. 

The  problem  is  to  optimize  the  firing  rate  (i.e.,  heat  transfer  rate)  for 
a  given  tube  weight. 

The  addition  of  fins  to  the  surface  of  a  tube  seriously  complicates 
the  calculation  of  heat  transfer.  Conduction  of  energy  within  the  fin  must 
be  calculated  simultaneously  with  the  radiative  and  convective  losses  at 
the  fin  surface.  Mitual  irradiation  between  adjacent  fins  and  between  the 

m 

fins  and  the  tube  surface  must  also  be  considered.  By  contrast,  the  radia¬ 
tive  and  convective  heat  transfer  rates  from  a  smooth  cylinder  are  independent 
of  each  other  and  may  be  evaluated  from  relatively  simple  expressions. 
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Finning  also  complicates  the  prediction  of  tube  strength  although 
it  should  be  noted  that  annular  fins  contribute  significantly  to  the 
load  carrying  capacity  of  the  tube.  For  a  given  ir.temal  pressure,  a 
properly  designed  finned  tube  has  a  smaller  wall-ratio  than  a  structurally 
equivalent  smooth  tube,  although  the  overall  tube  weight  per  unit  length 
is  slightly  higher  in  the  former  case.  Thus,  the  weight  penalty  incurred 
with  the  addition  of  annular  fins  is  less  than  the  actual  weight  of  material 
making  up  the  fins. 

The  problem  facing  the  designer  can  be  stated  as  follows:  for  a  given 
weight  addition  of  material  per  init  tube  length,  maximite  the  heat  dissi¬ 
pating  ability  of  the  outer  surface  while  simultaneously  maintaining  the 
structural  integrity  of  the  tube  under  the  specified  loading  conditions. 

The  heat  transfer  rate  and  the  tube  section  strength  must  be  considered  to¬ 
gether  since  they  are  related  through  the  maxinun  allowable  material  tempera¬ 
ture  and  the  tube  geometry  —  that  is,  the  tube  wall-ratio  and  the  spacing, 
height  and  profile  of  the  fins. 

A  desirable  design  tool  would  appear  to  be  a  computer  program  which 
accepts  as  input  such  basic  parameters  as  weapon  caliber,  internal  peak 
pressure,  material  properties  —  that  is,  the  thermal  conductivity,  surface 
radiative  emissivity  and  yield  stress,  all  at  the  maxinun  tube  temperature  — 
and  the  maxinun  tube  weight  per  unit  length.  Output  would  be  the  tube 
geometry  and  the  optimized  heat  transfer  rate  per  unit  tube  length. 

Information  in  two  areas  is  required.  First,  an  empirical  relationship 
mist  be  found  which  relates  tube  strength  and  material  yield  stress  to  an 
arbitrary  fin- tube  geometry.  This  information  can  be  obtained,  for  example, 
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through  a  parametric  photoelastic  study  of  various  fin-tube  sections. 
Secondly,  the  ability  to  predict  the  heat  transfer  rate  from  an  arbitrary 
tube  geometry  is  necessary.  Bnpirical  information  on  the  convective  heat 
transfer  coefficient  as  well  as  a  mathematical  model  of  the  combined  radia¬ 
tion-convection  problem  must  be  obtained.  The  latter  problem,  along  with 
a  cursory  examination  of  the  optimization  problem,  will  be  treated  in  this 
report.  The  photoelastic  and  convection  studies  will  be  the  subject  of 
future  proposed  work. 

PROBLEM  STATEMENT 

The  quasi-steady  nature  of  the  temperature  distribution  on  the  outer 
barrel  surface  permits  a  steady- state  heat  transfer  analysis  to  be  made. 
Rirthermore,  the  spatial  temperature  variation  in  the  axial  direction  can 
be  taken  into  account,  when  desired,  fay  subdividing  the  tube  into  sections 
along  its  axis  and  assigning  each  section  an  average  constant  temperature 
over  its  length.  The  segment  containing  the  temperature  maximum  is  of  pri¬ 
mary  interest  here. 

The  tube  geometry  can  be  approximated  as  shown  in  Figure  1.  An  actual 
mortar  design,  of  course,  would  have  fillets  at  the  intersection  of  the  fins 
and  the  tube  and,  although  this  would  be  an  important  consideration  in  regard 
to  tube  strength,  it  is  not  expected  to  substantially  alter  the  heat  transfer 
prediction.  Thus,  the  problem  reduces  to  an  analysis  of  the  combined  steady 
convective  and  radiative  heat  transfer  from  an  ensemble  of  annular  fins  of 
trapezoidal  profile  on  a  circular  isothermal  tube. 
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Convective  transfer  from  annular  fins  of  various  profiles  has  been 
summarized  by  Kraus  [1].*  Radiative  transfer  from  fins  of  uniform  profile 
was  studied  by  Sparrow,  Miller  and  Jons  son  [2] .  Their  analysis  included 
the  conplicated  radiative  interaction  between  a  fin  and  the  tube  base  as 
well  as  mutual  irradiation  between  adjacent  fins.  Edwards  and  (haddock  [3] 
extended  the  treatment  to  include  convection  but  limited  their  calculations 
to  specific  fin  materials  and  tube  teng>eratures.  Dent  [4]  treated  the  same 
problem  but  confined  his  work  to  a  method  of  solution  rather  than  a  presents* 
tion  of  design  information.  Recently,  Keller  and  Holdredge  [5]  presented 
results  for  radiative  transfer  from  a  single  annular  fin  of  trapezoidal 
profile.  Convection  and  incident  radiation  upon  the  fin  surface  were  not 
included  in  their  analysis.  The  present  study,  therefore,  represents  an 
extension  of  these  analyses. 

ANALYSIS 

Basic  Assumptions.  The  configuration  to  be  analysed  is  shown  in  Figure  2. 

The  nomenclature  is  shown  in  the  upper  portion  of  the  figure.  The  non* 
dimensionalized  variables  shown  in  the  lower  portion  of  the  figure  will  be 
defined  later. 

Because  an  ensemble  of  fins  on  an  isothermal  surface  is  assumed  only 
the  center-line  to  center* line  distance  between  any  two  adjacent  fins  has 
to  be  considered.  It  is  further  as  suited  that: 

1.  Conduction  of  heat  within  the  fin  is  one-dimensional  and  steady. 

2,  The  fin  and  tube  surfaces  are  gray  diffuse  reflectors  and  emitters 

*  Ntnbers  in  brackets  designate  references  at  end  of  report. 
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with  constant  emissivity,  c. 

3.  A  constant  convective  film  coefficient,  h,  prevails  over  all 
surfaces. 


4.  The  radiation  from  the  environment  can  be  characterized  as  coming 
from  a  black  body  source  at  temperature  Te*  equal  to  that  of  the  surround* 
ing  medium. 

5.  The  fin  base  and  tube  surface  have  the  same  constant  temperature, 


V 


6.  The  thermal  conductivity  of  the  fin  material,  km,  is  a  constant. 
Energy  Equation.  Consider  the  incremental  volume  shown  in  Figure  2  located 
between  r^  and  +  dr^.  Ch  the  fin  surface  this  volune  element  has  a 

A 

surface  area  dA^  -  r^dedr^/cos^  where  rjde  is  the  width  of  the  element  and 
dr^/cosat  is  the  height  (see  Figure  3).  Energy  is  conducted  into  the  element' 


at  r^.  It  loses  energy  by  conduction  at  rj  ♦  drj  and  by  convection  and 

A 

radiation  at  its  surface,  dAj.  An  energy  balance  gives 


d_ 

drl 


Vll6  ♦  (R0-rl)tm^l<Irl(rl)l  dedrl 


-  qjCr^rjdedrj  ♦  hlTjCrp-T^rjdedrj  (1) 

cos  #<  cos  •< 

The  term  on  the  left  represents  the  net  energy  conducted  into  the  element  per 
unit  time.  The  variable  fin  thickness  is  [6  ♦  (F0-ri)tan«<]  and  1%  is  the 


material  thermal  conductivity. 


The  first  term  on  the  right  represents  the  net  energy  lost  by  radiation 
from  the  element  per  unit  time.  qiC?])  is  the  radiative  heat  transfer  per 
unit  area  per  unit  time  at  r^«  The  second  tern  represents  the  convective 
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loss  due  to  the  temperature  difference  [Ti(ri)  -Te]  between  the  fin  and 
the  surroinding  nediun. 

Equation  (1)  is  a  second-order  ordinary  differential  equation  involv-  . 
ing  two  unknown  finctions,  T^(ri)  and  qi(r^).  An  appropriate  equation 
will  be  derived  in  the  next  section  for  the  radiative  flux,  qi(rj).  Of 
the  two  boundary  conditions  required  on  Tj(rj)  the  first  is  given  by 
assumption  (S)  above  which  states  that  the  fin  base  temperature  is  equal 
to  the  tube  temperature,  Tb,  or 

Tl(rl)  ■  Tb  at  rx  -  R£  (2) 

The  second  boundary  condition  is  obtained  by  equating  the  energy  conducted 
to  the  fin  tip  area,  6Pjde,  to  the  energy  lost  by  convection  and  radiation 
from  this  surface,  or  , 

-VTl(n)  *  h[Ti(ri)  -  Te]  fiRode 

drl  4  4 

♦  cfoTjCrj)  -oTe]  6RQde 

at  rx  -  1^  (3) 

A  similar  procedure  could  be  followed  for  the  temperature  distribution 
T2(r2)  in  the  opposite  fin,  but,  due  to  symmetry,  T2(r2)  is  identical 
to  the  function  T^r^)  when  r2  ■  r^,  hence,  the  resulting  equation  would  be 
redundant. 

Radiation  Equations.  Attention  will  now  be  directed  to  the  determination  of 

the  radiative  flux  term,  qi(r^),  in  equation  (1).  The  net  radiation  per  unit 

* 

area  per  uvlt  time  leaving  the  surface  ele.nent  dA^  is  equal  the  energy  emitted 
from  this  surface,  by  virtue  of  its  temperature,  less  the  energy  absorbed. 

The  latter  quantity  includes  contributions  from  the  surroundings,  the  opposite 
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fin  and  the  tube  surface.  It  also  includes  energy  originally  emitted 

* 

by  dAj  and  arriving  back  at  this  surface  after  reflections  off  the 
opposite  fin  and  tube  surfaces  when  these  surfaces  are  not  black  (i.e. , 
c<l).  The  complex  radiation  balance  is  most  easily  made  by  employing 
enclosure  theory  and  the  concept  of  radiosity  [see  ref.  6,  Chapter  3]. 

The  enclosure  to  be  analyzed  is  shown  in  Figures  2  and  3.  It  consists 
of  the  conical  surfaces  of  two  adjacent  fins,  designated  hy  the  subscripts 
"l”  and  "2",  the  cylindrical  base  surface  between  the  fins,  designated 
by  the  subscript  "b",  and  the  inner  surface  of  the  fictitious  cylindrical 
sheath  stretched  across  the  opening  shown  in  Figure  2  and  designated  by 
the  subscript  "e".  As  noted  earlier,  this  sheath  is  treated  as  a  black 
body  with  an  effective  temperature,  Te. 

The  term  "radiosity"  refers  to  the  total  energy  per  unit  area  per  unit 
time  leaving  a  surface.  It  is  equal  to  the  sun  of  the  diffusely  emitted  and 
diffusely  reflected  radiation  and  may  be  written  as 

Bitti)  ■  c^Cli)  ♦  (1-tjJHjdj)  (4) 

The  subscript  "i"  designates  the  i-th  surface  of  the  enclosure  and  the  vector 
r^  designates  the  position  of  an  element  on  this  surface  with  respect  to  some 
convenient  coordinate  system.  The  first  term  on  the  right-hand  side  of  eo. 
(4)  represents  the  rate  at  which  energy  is  emitted  per  unit  area  per  unit 
time;  is  the  emissivity  of  the  i-th  surface.  The  symbol  is  called 

the  "irradiation"  and  refers  to  all  radiation  incident  at  Ji  per  unit  area 
per  unit  time.  The  term  (1-ej)  represents  the  reflectivity  of  a  gray  surface 
and  the  product  (l-ej)Hi(ri)  is  the  energy  per  unit  area  per  unit  time 
diffusely  reflected  from  the  surface  at  r^. 


Since  the  radiosity,  CXi)  •  is  assumed  to  be  diffusely  distributed, 

the  concept  of  the  angle  factor  can  be  employed  [see  ref.  6,  Chapter  4]. 

In  particular,  the  notation  dF(r^,  Xj^dA^-dAj  b*  use<*  t0  <iesiSnate 

the  fraction  of  the  energy  leaving  surface  element  dA?  located  on  the  i-th 

surface  at  r^  and  arriving  at  dAj  located  on  the  j  -  th  surface  at  jj  *  ^h® 

asterisk  on  dA^  is  used  to  denote  a  "second-order"  differential  area  such 
*  * 

as  dAj  in  Figure  3;  that  is,  dA^  ■  r^dedr^/cos  •(.  When  the  asterisk  is 
omitted,  a  "first-order"  differential  area  is  implied;  for  example, 
dA2  •  2irr£  dr 2 /cos  «( . 

The  notation  F(i i)dA*-A-  wil1  **  use^  to  designate  the  fractior  of  the 

^  * 

energy  leaving  surface  element  dA^  which  is  intercepted  by  the  entire  surface 
Aj.  The  following  "reciprocity  rules"  are  derived  in  ref.  6  and  will  be  used 
in  the  development  below. 

dAiF^-i^d^-Aj  "  AjF(Ji)Aj-dAj  (6) 

The  expression  for  the  energy  incident  on  dAj  will  now  be  derived. 
Consider  first  the  contribution  from  the  opposite  fin.  The  total  energy  per 
unit  time  leaving  the  ring  element  dA2  located  at  r2  in  Figure  3  is  B2(r2)dA2. 
Of  this  amount  the  fraction  dF(r_,  r,),.  * 

rPdA2-dAj 


arrives  at  dAj  or 


The  total  irradiation  of  dA^  from  the  opposite  fin  is  found  by  summing  the 
contributions  of  all  of  the  ring  elements  on  A2  with  the  result 
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(  B2  Cr2>dA2dF(x‘2  *  rl^dA2-dA* 

h 


Specializing  eq.  (5)  to  this  situation  gives 

dAgdFO^,  -  dAjdFCrj,  r2*c!A*-dA2 

Using  this  result  in  the  above  integral  gives 

A 

B2^r2)dF^rl»  r2^dA*-dA2 

A 

where  dA*  has  been  removed  from  under  the  integral  sign  since  it  is  indepen¬ 
dent  of  the  integration  variable  r2. 

A  similar  procedure  can  be  followed  to  calculate  the  irradiation  of 

* 

dAj  by  all  ring  elements  dA^  on  the  tube  surface  with  the  result 
L  -  (R^R^tan^ 
dAj  (^(z^dFCrj,  *1,)^*.^ 


(R0  -  tytan* 

B^(z^)  is  the  radios ity  of  the  ring  element  d\  located  at  zj,  on  the  tube 
surface.  Note  the  asymmetric  origin  of  the  z^ -coordinate  in  Figure  2  end 
note  also  that  the  tube  surface  extends  from  -  (RQ-R^tan^  to  z^,  ■  L  - 
(Ro-Ri)tan^  . 

Since  the  surface  A^  is  assumed  to  behave  as  a  black  body  with  uniform 
temperature,  Te,  its  emissivity  is  unity  and,  by  eq.  (4),  its  radiosity  is 


4 

simply  Be  ■  oTe,  independent  of  position.  The  total  energy  emitted  by  Ag 
is  oT^Ag  and  the  fraction 

is  intercepted  by  dA^.  Using  the  reciprocity  rule  in  eq.  (6)  this  quantity 
can  also  be  written  as 

Since  surface  "1"  is  convex  it  does  not  irradiate  itself,  so  the  total 
energy  per  unit  time  incident  on  dA*  is  the  sun  of  the  three  contributions 
above.  Dividing  this  sun  by  dA^  the  irradiation  per  unit  area  per  wit 
time  is  obtained,  that  is,  H^(r^).  Thus, 

Hl(rl)  -  \  B2(r2)dF(r1,  ^dA^-dAj 

/ri 

L  -  (V  RiK*n*  (?) 

♦  *  <,T2F<rl>dAj-AB 

'■Ro"R1)tan«( 

The  net  local  radiation  flux  is  given  by 
4 

qiCli)  “  H<*Ti(ri)  -  ciHi(ji)  (8) 

The  second  tern  on  the  right-hand  side  of  eq.  (8)  represents  the  fraction 
of  the  incident  energy  which  is  absorbed;  by  the  gray  body  assumption  the 
absorbtivity  is  equal  to  the  emissivity,  e^. 

Eliminating  H^(n)  in  eqs.  (4)  and  (8)  gives 
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(9) 


ciBi(Xi)  ■  cioTjdp  -  (1  -  e^q*^) 

Using  this  expression  to  eliminate  B2(r2)  and  B^Cz^)  in  eq.  (7)  and  sub¬ 
stituting  the  result  into  eq.  (8)  gives 

»  n 


qi Crx)  ■  coTj(rj)  -  \  [eoT2(r2)  -  (l-c^f1^)  JdF(rl»  ^dA^-dA- 

h 

'L  -(V^*®1* 


-IKo-KiJzan^ 

-  ([c.T«  -  (l-eJ^Cu.ndFCr!,  «„),**„.  -  ^F(rl)dA*.v 

'(R^tan^ 

(10) 

Due  to  symmetry  the  expression  for  q2(r2)  can  be  obtained  by  simply 
interchanging  the  subscripts  "1"  and  "2”  in  the  above  equation.  Also,  note 
that  the  tube  surface  temperature,  Tjj,  is  not  a  function  of  position  (see 
assumption  (5)  above). 

An  analogous  procedure  yields  the  following  expression  for  the  base 
radiative  flux: 

A> 


%(*b)  -  «1b  '  U'oTta)  -  (l-OqitrjJdFdb.rj)^.^ 


A 

(  [c«T*(r2) 
\ 


{l-c)q2(r2)]dFCxb. 


-  eoT, 


(11) 


Note  that  while  the  tube  emits  energy  uniformly  over  its  surface  by  virtue 
of  its  constant  temperature,  the  irradiation,  denoted  by  the  last  three  terms 
in  eq.  (IP,  is  a  function  of  position.  Finally,  because  the  tube  is  a  convex 
surface  it  does  not  irradiate  itself. 


Eqs.  (10)  and  (11)  can  be  simplified  considerably  by  taking  advantage 
of  certain  symmetries  in  the  problem  and  through  the  use  of  angle  factor 
algebra.  As  noted  earlier,  due  to  symmetry,  the  functions  T2(rg)  and 
q2(rz)  are  identical  to  the  functions  T^(rj)  and  q^(r^)  when  r2  ■  r^.  Also, 
the  angle  factor  dF(ib,  ^dA^-dAj  "“V  **  replaced  by  dFtL-z^rj)^.^ 
when  r2“  r^.  The  latter  angle  factor  is  evaluated  at  (L-i^)  rather  than 
z^  because  the  fins  are  not  symmetrically  located  about  the  z^  origin  (see 
Figure  2).  Finally,  the  following  angle  factor  relationships  apply: 

F<rl>dA^  ♦  (  ^l^djydAj  4  (  "fria^-dAt  "  1  ™ 

(R0-Ri)tan  k 

R 0  ^ 

F(tb)dA^-Ae  *  (  "(‘b.nJdAj-dAj  4j  "('b.  r2)dA^-dA2  *  1  (13> 

Ri 

The  three  terns  an  the  left-hand  side  of  eq.  (12)  represent,  respectively, 
the  fraction  of  the  energy  leaving  dA^  which  is  intercepted  by  the  sheath 
area,  A^,  the  area  of  the  opposite  fin,  A2,  and  the  tube  base  area, 

Since  these  three  surfaces  and  Aj,  which  does  not  see  itself,  foim  an 
enclosure,  the  above  three  fractions  must  sun  to  unity  as  indicated  in  eq. 
(12).  A  similar  interpretation  applies  to  eq.  (13). 

Adopting  the  notation 

*12 (rl*  r2>  -  dF<rl»  r2>dAi-dA2 

Klb(*l.  lb)  *  dF(rl>  ^dAj-dAj,  /dlb 
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KfalCV  rl>  “  ^JdA^-dAj  /drl 

and  using  the  above  results,  eqs.  (10)  and  (11)  become 

R° 

<ll(ri)  -  etoTjtr!)  -oT,]  -  (  U[ot}(0  -°T*]  -(l-cJq^O)  r12(r1(t)d{ 

'Ti  . 


Ri 

L  -  (R^R^tan* 

-  (  <e[oT|J  -  oT*]  -  (l-ejqb(0>  K^CTj,  C)dC 

XR0-Ri)tanK 

% 

%(*!,)  -  e[oT||-oT*]  -){c[ot}(0  -oT*]  -  (l-Oq^)) 

Ri 
Ro 


(14) 


i 


(c(oTj(t)  -otJ]  -  (l-c)q1(0)  ^(L-^.Odt 


(15) 


As  is  typical  in  problems  of  this  type  the  kernel  K^O^fO  in  the 
second  integral  of  eq.  (14)  becomes  indeterminate  at  the  upper  limit 
[L  -  (R0-R^)tan«i]  when  r^  ■  R^,  that  is,  at  the  intersection  of  the  fin 
and  the  base  cylinder.  The  difficulty  in  evaluating  the  integral  numerically 
at  this  point  can  be  circumvented,  however,  by  following  a  procedure  given 
by  Donovan  [8]  for  analytically  determining  the  integral  in  the  limit  as 
T}  -*■  R^.  The  details  are  too  lengthy  to  be  reproduced  here  and  only  the 
result  will  be  quoted.  Thus 

limit 
rl  ♦  % 


L  "  CVRi)1®1^ 

{e[oTj-oTj]  -(l-c)qb(C)}Klb(r1,OdC  - 


(Ro-Ri)tanK 


-  1c[«tJ|-oT*]  -  (l-Oq^l  -  (R0-Ri)t«n«OHl-siM)/2 

A  similar  difficulty  occurs  at  the  lower  limit  in  the  first  integral 
in  eq.  (IS)  as  z^  ■*>  L  -  (R^R^tanof.  Analytical  evaluation  gives 

fr 

limit  \U[oT}(r.)-0T*]  -(l-€)q,U))ICblCtb,Oit 

-L-(R0-Ri)tanon 

-  <c[oT|J-<,T*]  -  (l-c)q1(Ri)Kl-sin.O/2 

4  4 

where  eq.  (2)  has  been  used  to  replace  T^(R^)  by  in  the  first  term  on 
the  right-hand  side  of  this  equation. 

The  second  integral  in  eq.  (15)  can  be  handled  in  the  same  manner  as 
>b  -*•  (Ro-Ri)tan*  or  advantage  can  be  taken  of  the  symnetry  of  %(zb)  about 
the  point  ■  L/2  by  noting  that 

<lbC(Ro-Ri)tan-0  ■  qb(L  -  (Ro-R^tan^) 

Che  final  simplification  resulted  after  evaluation  of  the  angle  factors 
revealed  that 

Klb(rl»  rb)  *  (Ricos  H/r^KbjCzb,  rx) 

Expressions  for  Kj2(rl»  r2)  ^  ^blC^*  rl)  arc  giv®1  in  Appendix  A. 

Non-dimensionalized  Equations.  The  following  dimensionless  variables  were 
used  to  non-dimensionalize  the  governing  equation. 

Cl  "  n  »  ?2  "  £2  »  Tb  "  fb 

%  Ro  Ro 

NL  -  _L 

Ro 


nr-  % 

Ro 


I 
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» 


6*  -  6 

R^ 


TiCCi)  -  TiCrj)  , 


Te  "Te 


^lCCi)  "  <U(ri) 


Qb(Tb)  "  %C*b) 


ot: 


oT, 


T 


N. 


h*  - 


n2  t3 

Ri  0  Tb 


cos*(a  ♦  (Ro-R^tanA) 


hRi 


V„,cosA(«  ♦  (R^R^tan*) 

With  this  notation  the  final  form  of  the  radiation  equations  becomes: 


Qi(Nr)  -  t[l-T*4]  -  \  (t[T*4(?)  -  T*4]  -  Cl-t)Qi(0)K12(NR>  C)dt 


% 


*4 


-  (t[l-Te']  -  (l-c)Q^(NL-Cl-Np) tan* ) ) {1-sin ^ )/2 


(16a) 


W 


-  |<«rrJ4(0-T*4)  -  (l-OQjCOUtnCfi.Odt 


Nl  -  (l-NR)tan*( 

-  NrcosK  \(c(l-T*4]  -  (l-eJl^CO^jU,  ft)d?  (16b) 

/(l-NR)tan< 

Qb(N,-(l-NR)tan.O  -  e[l-T*4]-U[l-T*4J-(l-e)Q1(NR))(l-sin.O/2 

-((e(T*4(0-T‘4]  -  (l-c)Q1(0)*bl((l*NR)tan.U)dC 


% 


(17a) 
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(17b) 


<^((l-NR)t*n<)  -  <^)(Nl-(l-NR)t«n.O 

VTb>  ■ 


£• 


{c[Tl4(0-Te4]  ■  (l-e)Qi(0>[Kbl(Tb.  0  ♦  ^10^.0  ]<U 

(17c) 


The  conduction  equation  and  its  boundary  conditions  become 


d 

d(l 

where 


P(Ci)  ^lCCl) 
Oft 


ft  (NcQl(ft)  *  hMTiCft)  -Ten 


(18a) 


N) 


P(ft)  *  fj«»  ♦  (l-fl)MMt) 

6*  ♦  (1 -Nr)  tan/ 

T*(fl)  -1  at  f!  -  Nr  (18b) 


C  ■  -  dT^(f^) 

dfl 

-  (h*IT*(ri)  -  T*]  *cN.lT*4(fl)-T*4]) 

x  [«*  ♦  (l-Hj^tan*  ]cos<  at  1  (18c) 

>5 


Method  of  Solution.  Following  a  technique  suggested  by  Glauz  [7]  and 
successfully  applied  by  Donovan  [8]  to  the  related  problem  of  inifoim  fins 
on  a  plane  wall,  eq.  (18a)  was  transformed  into  an  equivalent  integral  equa* 
ticn  by  direct  integration  and  application  of  the  boundary  conditions.  The 
result  is 
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Tl«y 


_ 1_ 

PtCi) 


X1  A 

(Nr-^)  (  ♦(Ode  ♦  f#(0<U  ♦  nr 


4 

(nr-(°i)4,c 


6*  ♦ 


(l-N^tan  •< 


where 


♦CO 


(C/Nr)  (WO*  hVl(O-Tel) 


(19) 


♦(O  -  *  (l-2C)t*n*(l  *(Nr-0*(E) 

[«’  ♦  (3 -Nr)  tan 

The  constant  C  is  given  in  eq.  (18c)  with  ^  »  1. 

The  advantages  of  the  integral  formulation  in  eq.  (19)  over  the  differ¬ 
ential  equation  in  eq,  (18a)  are:  (1)  the  bomdary  conditions  are  included 
in  the  equation  and  (2)  the  equation  is  in  the  same  form  as  the  radiation 
equations,  a  form  which  readily  lends  itself  to  an  iterative  type  solution. 

fyecific  details  regarding  the  solution  of  eqs.  (16) ,  (17)  and  (19)  for  the 

• 

distributions  of  Qi(fi),  Qjj(t^)  and  Ti(fj)  are  given  in  Appendix  B. 

Total  Heat  Transfer  and  Fin  System  Effectiveness.  The  heat  transfer  rate 
from  the  finned  tube  section  between  the  two  vertical  center-lines  in  Figure 
2  is  equal  to  the  sun  of  the  radiative  and  convective  contributions.  The 
radiative  contribution,  R,  is 


j/%  L  "  (Ro-Ri)tan< 

qR  *  2  \  ql(rl)  2,rldrl  *  (Vv  2,Ridjb 

/rj  c05*< 

♦  2  (2»R0«)  [taff^CRo)  -  Te4)] 


(20) 


The  integrand  under  the  first  integral  in  eq.  (20)  is  the  product  of  the 
local  radiative  flux,  qi(r^),  aid  the  area  of  the  ring  eleamnt,  dA},  in 
Figure  3.  Integrating  over  the  entire  fin  surface  area  and  Multiplying 
the  result  by  the  factor  2  to  take  into  account  the  contribution  the 
apposite  fin  surface  ^  gives  the  net  radiative  loss  from  the  two  conical 
surfaces.  The  second  integral  gives  the  net  radiative  loss  from  the  base 
surface  between  the  fins.  The  last  tens  represents  the  radiative  loss 
from  the  two  fin  tips  of  area  (2*1^6)  at  tenywrature  T^CRq). 

The  convective  loss,  is  written  similarly  as 

%  “  2  \  hCTjCrj)  -Te)  2wr1dr1 
*4.  cos  *< 

♦  (2wRi)(L-2(R0-Ri)tan<)h(Tb-Te)  ♦  2(2»R0«)h(T1(R0) -T^  (n) 

where  the  second  term  represents  the  convective  loss  from  the  isothermal 
tube  surface. 

A  rn— nn  method  of  assessing  the  enhancement  in  total  heat  transfer 
rate,  q^,  due  to  the  addition  of  fins  is  to  compare  this  value  with  the 
rate  of  heat  transfer  from  the  uifinned  cylinder,  qQ.  Hie  latter  quantity 
is  given  by 

q„  •  2.R1(L-2«)[co(t5-tJ)  .  h(Tb-T,)J  (22) 

Hie  comparative  index  is  called  the  "fin  system  effectiveness"  by  Donovan 
[8]  and  is  written  as 


1 


Fins  are  effective  in  increasing  the  heat  transfer  rate  when  the  fin 
configuration  md  thermal  properties  are  such  as  to  produce  values  of  3 
greater  than  inity,  Fro*  the  nmufacturing  point  of  view,  3  should  be 
substantially  above  mity.  Systens  with  3  less  than  inity  must  be  avoided 
because,  in  this  case,  the  fins  have  an  insulating  effect  on  the  base  tube. 
In  dimensionless  form  the  fin  system  effectiveness  is  given  by 

5  -  Nc  q,  ♦  h%,  (24) 

% 

’*"*  A  /Nl  -  (l-ftytm* 

Or-  2  (  <jl(fi)(jd<i  ♦  N„  (  (^(Tb)d,b  ♦  2S\[tJ4(1)-T*41 

\  «*<  '(l-N^Un* 

Qc  *  2  ♦  NR(l-VtNL-2q-NR)t»l^1  ♦  2«Vl(l)-T.l 

Nn  COS  •< 


Qo  -  Wl  ♦  26  )[cNc(l-Te4)  ♦  h  (l-Te)J 

The  computer  program  listed  in  Appendix  C  calculates  the  distributions 

* 

of  ^1^1^»  W  T^irou8^  the  above  equations  it  then  computes  3, 

With  a  little  algebra  it  can  be  shown  that  the  actual  heat  transfer  rate 
is  related  to  3  and  by  the  expression 

<*!•  *  (ZwjQo/NRjCRoTbl^cosOCa*  ♦  (l-N^tan*)  (25) 


APPLICATION  OF  THEORY  TO  TOE  EESIGN  OF  A  MORTAR  TUBE 
Input  Data.  As  an  application  of  the  theory  the  heat  transfer  character¬ 
istics  of  the  fin  configuration  used  on  the  XM144E1  mortar  tube  (see  Figure 
1)  were  calculated.  The  ambient  air  temperature ,  Tft,  was  taken  as  100*F 
in  all  of  the  calculations.  The  emissivity,  c,  the  metal  thermal  conductivity 
and  the  tube  temperature,  Tb,  were  varied  as  shown  in  Table  I  below. 

Table  I 
c  *  0.75,  1.0 

-  10,  25  Btu/hr-ft-'F 
Tb  -  800*F  -  1200#F 

The  emissivity  values  are  typical  of  the  coatings  used  on  mortars  operating 


in  this  tenperature  range.  The  low  value  of  is  typical  of  a  high  nickel 
content  alloy  while  the  higher  value  corresponds  to  a  higher  iron  content. 

The  three  parameters  —  c,  md  Tb  —  are  characteristics  of  a 
particular  material;  the  maximum  allowable  tube  tenperature  is  related  to 
the  mechanical  properties  of  the  tube.  It  is  clear  before  any  results  are 
presented  that  the  material  with  the  highest  values  of  this  set  of  parameters 
will  give  the  highest  heat  transfer  rate  —  the  present  calculations  are 
intended  to  show  the  relative  importance  of  each  of  the  parameters.  The 
effect  of  changes  in  geometry  for  given  e,  k„,  and  *ill  be  considered  in 
a  later  section. 

Calculation  of  the  Dimensionless  Convection  Coefficient,  h*.  The  dimension- 


where  the  product  (Lkft)  has  been  introduced  in  the  mine ra tor  and  denominator 
(kg  is  the  thermal  conductivity  of  the  air  and  L  is  defined  in  Figure  2) . 

Hie  dimensionless  grouping  (hL/kg)  is  recognized  as  the  Nusselt  number. 

Correlations  for  free  convection  heat  transfer  from  finned  tubes  are 
given  by  Knudsen  and  Pan  [9]  and  Nwizu  [10].  Both  correlations  are  re¬ 
stricted  to  straight  «■  0)t  isothermal  (T,(r,)  ■  T*)  fins  on  a  horizontal 
cylinder.  In  addition,  their  tests  were  limited  to  tubes  with  diameters  of 
1-1/2"  or  less  and  a  maximum  temperature  of  3S0*F.  Although  the  short  fins 
shown  in  Figure  1  do  operate  in  a  nearly  isothermal  condition  the  limitations 
an  surface  temperature  and  particularly  tube  orientation  raise  considerable 
doubt  as  to  the  applicability  of  their  work  to  mortar  operating  conditions. 
The  geometric  limitations  also  introduce  some  uncertainty.  Thus,  while  these 
correlations  should  not  be  applied  in  an  actual  mortar  design  situation  they 
do  serve  the  purpose  here  of  relating  the  convection  coefficient  to  the 
taupe rature  dependent  air  properties  and  to  the  fin  geometry.  The  need  for 
further  experimental  work  in  this  area  is  evident. 

Over  the  range  of  variables  of  interest  here  Nwizu' s  correlation  can 


be  written  as 


hL  ■  0.058 


gat2  l3  cvr«) 

LX-2 


0.4S 


(26) 


where  0  is  the  coefficient  of  thermal  expansion  (1/*F) ,  f  is  the  density 


(itfo/ft3),  g  is  the  acceleration  due  to  gravity  (32.2  ft/sec2)  and  /*i  is 
the  viscosity  (lb^/ft-sec.).  All  air  properties  are  evaluated  at  the 
temperature  Ob  *  Tc)/2. 
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Heat  Transfer  Results,  The  results  of  the  above  calculation  are  presented 
in  Figures  4  and  S.  Of  the  two  sets  of  curves  plotted  in  Figure  4  the 
upper  set  represents  the  total  heat  transfer  per  wit  tube  length;  that  is, 

“  %+  <*c  (27) 

S  S 

where  and  are,  respectfully,  the  radiative  and  convective  portions  of 
the  total  heat  transfer  (see  eqs.  (20)  and  (21))  from  the  fin- tube  section 
of  axial  length  S  -  L  ♦  26  (see  Figure  2).  The  lower  set  of  curvv  in 
Figure  4  represents  the  radiative  contribution  a?  one  where  qR  ■  qR/S. 

Of  the  three  parameters  shown  in  the  figure  the  base  temperature,  Tb, 
has  the  dominant  effect  on  the  heat  transfer.  This  is  primarily  due  to 
the  tremendous  increase  in  radiative  transfer  with  temperature.  Consider 
the  uppermost  curve  in  each  set,  that  is,  c  -  1.0  and  -  25  Btu/hr-ft-  °F. 

At  Tb  “  800  "F  radiation  accounts  for  roughly  561  of  the  total  heat  transfer. 
Increasing  the  base  temperature  to  1200  "F  more  than  doubles  and,  at  this 
temperature,  radiation  accounts  for  over  701  of  the  total. 

These  results  can  be  explained  by  noting  that  the  potential  for  con¬ 
vective  transfer  is  (”%-  Te)  while  that  for  radiative  transfer  is  (l£- 
Increasing  Tb  from  800  IF  (1260  H)  to  1200  T*  (1660  *R)  while  keeping  T^  * 
100°F(560*R)  increases  the  tenperature  potentials  for  convective  and  radia¬ 
tive  transfer  by  roughly  578  and  2081,  respectfully.  The  actual  calculations 
show  a  538  increase  in  convection  and  a  2008  increase  in  radiation;  the 
discrepancy  is  due  to  the  tenperature  drop  along  the  fin  which  reduces  the 
average  potentials  for  heat  transfer  somewhat,  lie  magnitude  of  the  radiative 
contribution  indicated  in  Figure  4  justifies  the  tedious  calculations  required 
to  obtain  it. 
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Reducing  the  surface  emissivity  from  unity  to  c  ■  0.75,  a  2SI 
decrease,  reduces  the  radiative  transfer  by  only  about  half  this  amount. 

Ids  surprising  result  is  explained  by  the  so  called  "cavity  effect". 

When  e  is  unity  the  only  radiation  which  escapes  from  the  fin-tube  cavity 
is  that  which  is  emitted  by  a  warm  surface  in  the  direction  of  the  fictitious 
cold  surface  Me"  in  Figure  2.  All  other  radiation  is  completely  absorbed 
by  the  remaining  black  surfaces  of  the  enclosure.  Decreasing  the  emissivity 
decreases  the  emissive  power  of  each  elemental  area  in  the  enclosure  but  at 
the  same  time  it  increases  the  probability  of  escape  of  radiation  from  the 
enclosure  by  the  process  of  multiple  reflections.  Thus,  radiation  emitted 
by  one  warm  element  in  the  direction  of  a  second  warm  element  may  be  reflected 
by  the  latter  surface  in  the  direction  of  surface  "e"  and  escape  from  the 
system.  The  net  result  is  a  decrease  in  net  radiative  loss  from  the  enclos¬ 
ure  which  is  substantially  less  than  the  decrease  in  emissivity.  The  converse, 
of  course,  is  also  true  —  increasing  the  emissivity  does  not  result  in  a 
proportionate  increase  in  radiative  transfer.  Therefore,  if  a  particular 
coating  is  found  to  be  satisfactory  from  the  viewpoint  of  durability  and  has 
an  emissivity  of  0.8-0, 9  at  the  desired  operating  temperature ,  very  little 
will  be  gained  by  searching  for  a  coating  with  a  higher  emissivity. 

The  thermal  conductivity  of  the  fin  material,  is  seen  to  have  a 
very  small  effect  on  the  overall  heat  dissipation  rate  in  Figure  4.  This 
is  far  from  the  general  state  of  affairs,  however,  as  will  be  seen  later 
when  the  effect  of  geometry  is  considered.  For  the  short  fins  under  con¬ 
sideration  here,  most  of  the  resistance  to  heat  transfer  is  encountered  at 
the  fin  surface.  The  temperature  drop  along  the  fin  amounts  to  only  a  few 
percent  for  all  of  the  cases  in  Figure  4. 
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Fin  System  Effectiveness,  The  fin  system  effectiveness,  3,  was  defined 
earlier  as  the  ratio  of  the  heat  transfer  rates  from  the  finned  and  smooth 
cylinders,  respectfully,  or 

3  •  It  ‘  2t  (28) 

«0  % 

A  plot  of  3  is  shown  in  Figure  5  for  the  fin  system  under  consideration. 

The  decrease  in  3  with  increasing  base  temperature  results  from  the  fact 
that  the  radiative  and  convective  contributions  to  increase  exactly  in 
proportion  to  the  potentials  and  (T^  -  Te),  respectfully,  whereas, 

as  noted  earlier  in  the  discussion  of  Figure  4,  the  average  temperature 
potentials  in  the  calculation  of  q-j  are  somewhat  lower  due  to  thermal  resis¬ 
tance  within  the  fin.  Thus,  the  ratio  in  eq.  (28)  is  a  decreasing  function 
of  temperature  and  the  fin  system  becomes  less  effective  at  higher  tempera¬ 
tures. 

Increasing  the  thermal  conductivity,  k^,  produces  a  relatively  small 
increase  in  3.  lie  increase  is  due  to  the  fact  that  energy  is  more  easily 
transported  to  the  outer  regions  of  the  fin  resulting  in  higher  average 
temperature  levels  and  a  more  effective  utilization  of  the  entire  fin.  In 
very  tall  fins  or  with  poor  conductors  the  thermal  conductivity  can  have  a 
much  more  pronounced  effect  on  j. 

The  surface  emissivity,  e,  has  a  strong  effect  on  3.  Decreasing  c  by 
2SI  decreases  the  radiative  contribution  of  the  smooth  cylinder  by  the  same 
amount  but,  due  to  the  cavity  effect  noted  earlier,  the  decrease  in  the 
finned  tube  heat  transfer  rate  is  only  about  half  this  amount.  Therefore, 
the  ratio  in  eq.  (28)  increases  with  decreasing  c.  Thus,  the  addition  of 
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fins  is  more  effective  in  augmenting  heat  transfer  in  systems  with  lew 
emissivity. 

With  regard  to  the  mortar  tube  configuration  shown  in  Figure  1  the 
results  in  Figure  S  indicate  that  the  addition  of  fins  represents  an 
augmentation  in  heat  transfer  on  the  order  of  501,  assuming  a  relatively 
high  emissivity.  This  increase,  while  significant,  still  leaves  a  wide 
gap  between  the  maximum  permissible  firing  rate  and  the  desired  rate;  that 
is,  the  rate  at  which  an  individual  can  load  the  weapon.  Assuning  a  material 
has  been  chosen,  that  is,  assuning  e,  and  have  been  specified,  the 

only  recourse  left  to  the  designer  is  to  alter  the  fin  geometry.  Ihis  will 
now  be  considered. 

Alternate  Fin  Configurations.  In  general,  the  designer  must  stay  within 
a  certain  weight  limitation  in  choosing  alternate  fin  configurations.  The 
weight,  or  equivalently,  the  volune  of  metal  per  unit  length  which  the  fins 
contribute  is 

V  -  f  4w 

Vs 

In  ?>e  present  set  of  calculations  V  was  kept  constant  and  equal  to  V0  which 
was  calculated  for  the  fin  configuration  shown  in  Figure  1;  that  is,  all  of 
the  alternate  designs  considered  here  have  the  same  weight  per  unit  length. 

In  addition,  the  fin  tip  thickness,  6,  and  the  tube  radius,  R^,  were  also 
fixed  at  the  values  shown  in  Figure  1.  The  fin  spacing  S  was  given  the 
values  Sq/2,  So.  3Sq/2  and  2S0  where  S0  is  the  fin  spacing  in  Figure  1. 

The  angle  •<  was  taken  as  0*,  2*,  4°,  6°  and  8°.  The  remaining  variable,  Rq» 
was  calculated  from  eq.  (29)  with  V  ■  V0.  The  resulting  twenty  configurations 
are  shown  in  Figures  6a  -  6d. 
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It  is  noted  in  the  figures  that  for  a  given  fin  spacing,  S,  decreasing 
•C  results  in  taller  fins  and  a  greater  surface  area  per  unit  tube  length. 
Increasing  the  fin  spacing  for  a  given  <(  also  increases  the  fin  height 
but  the  surface  area  per  unit  length  does  not  change  appreciably  because 
the  number  of  fins  per  unit  length  -decreases  with  increasing  S. 

Effect  of  Geometry  on  Heat  Transfer,  'he  total  heat  transfer  rate  per 
unit  length,  <5T»  was  evaluated  at  eight  sets  of  the  parameters  e,  and 
for  each  of  the  twenty  configurations  shown  in  Figures  6a-6d.  The 
results  are  presented  in  Figures  7a-7h.  Each  figure  gives  the  results  for 
one  set  of  e,  and 

Consider  Figure  7a.  Each  point  on  the  curved  surface  represents  the 
percentage  change  in  %  for  that  geometry  (S  and  *<  )  relative  to  qj  for 
the  XM144E1  fin  configuration  shown  in  Figure  1  at  the  values  of  c,  and 
Tb  listed  at  the  top  of  the  figure.  The  value  of  qT  given  in  the  figure 
title  corresponds  to  the  XM144E1  fin  configuration.  Thus,  the  fin  configura¬ 
tion  with  •  0*  and  S  -  2S.  (see  Figure  6d  for  this  geometry)  dissipates 
approximately  101  more  heat  than  the  XM144E1  mortar  or  <JT  »  1.1x0.870x10* 
Btu/hr-ft  ■  0.957x10*  Btu/hr-ft.  When  ■  8*  and  S  ■  S*/2  (see  Figure  6a) 
the  heat  transfer  rate  is  roughly  71  below  that  for  the  XM44E1  mortar  con¬ 
figuration.  Seventeen  of  the  twenty  configurations  in  Figures  6a-6d  are 
seen  to  dissipate  more  heat  than  the  configuration  shown  in  Figure  1. 

Comparing  Figures  7a  and  7b  shows  that  increasing  the  emissivity  causes 
a  general  flattening  of  the  surface  although  the  reduction  is  not  large. 

The  sanw  effect  is  noted  in  the  other  three  pairs  of  figures  when  the  emis¬ 
sivity  is  the  only  parameter  changed,  he  increase  in  c  increases  the  emis¬ 
sive  power  of  each  elemental  area  of  the  fin  surfaces  which,  in  turn, 
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increases  the  teqwrature  drop  along  the  fin  to  balance  the  energy  demand, 
fils  results  in  a  less  effective  utilization  of  the  outer  portions  of  the 
fin  and  the  effect  is  most  apparent  in  the  taller  fins  (#<  ■  0*).  Increas¬ 
ing  the  temperature  level  also  increases  the  radiative  transfer  and  has  the 
same  effect  on  the  shape  of  the  surface  as  c.  Compare  Figures  7a  and  7c 
or  Figures  7b  and  7d,  for  example. 

Increasing  the  thermal  conductivity,  1^,  markedly  inproves  the  heat 
transfer  because  of  the  more  efficient  transport  of  energy  to  the  outer 
portions  of  the  fins.  The  effect  is  strongest  as  the  fin  height  increases 
(decreasing  •{ ) .  Compare  Figures  7a  -  7d  with  Figures  7e  -  7h  respectively. 

When  the  fin  spacing  is  small  (S  ■  Sq/2)  ,  decreasing  «<  increases  the 
heat  transfer  because  of  the  increase  in  fin  surface  area.  At  large  fin 
spacings  (S  ■  ZSq)  the  heat  transfer  rate  at  first  increases  with  decreasing 
,  reaches  a  maximum  and  then  may  decrease  slightly,  as  in  Figures  7a  -  7d, 
The  decrease  can  ba  attributed  to  the  difficulty  in  transporting  energy  to 
the  outer  portions  of  the  taller  fins  («<  ■  0  ^  when  the  fin  thermal  con¬ 
ductivity  is  low. 

Concluding  Remarks.  Once  a  material  is  chosen,  that  is,  after  the  parameters 
c,  kjj,  and  lb  ftre  specified,  a  surface  such  as  those  in  Figures  7a  -  7h  can 
be  constructed.  The  problem  is  then  to  choose  the  highest  point  on  the 
surface  consistent  with  the  structural  integrity  of  the  tube.  The  importance 
of  being  able  to  predict  tube  strength  for  an  arbitrary  fin  configuration  is 
apparent. 

All  of  the  surfaces  seem  to  have  a  maxinun  in  the  interval  S0/?<S<S0; 
that  is,  the  analysis  seems  to  favor  small  fin  spacings.  This  conclusion 
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might  be  filtered  somewhat  if  the  orientation  of  the  mortar  tube  with 
respect  to  the  horizontal  could  be  taken  into  account  in  the  calculation 
of  the  convection  coefficient.  Wider  fin  spacing s  that  offered  less 
impedance  to  the  vertical  air  flow  would  probably  be  favored.  Since  the 
large  radiative  contribution  is  independent  of  the  tube  orientation  the 
results  are  probably  not  significantly  in  error,  however. 

Finally,  it  should  be  emphasized  that  all  of  the  heat  transfer  results 
are  based  on  one  tube  weight  per  isiit  length,  one  fin  tip  thickness,  6.  one 
tube  radius,  R^,  and  it  was  assuned  that  the  fin  material  had  a  conductivity 
on  the  order  of  that  for  steel  or  steel  alloys  —  as  opposed  to  a  better 
conductor  such  as  aluminum.  Thus,  this  study  should  be  viewed  more  as  a 
qualitative  analysis  or  method  of  interpretation  of  the  type  of  results 
which  can  be  obtained  with  the  computer  program  in  Appendix  C  rather  than 
a  quantitative  recommendation  for  future  mortar  designs. 
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(DN&USDNS  AND  REOfENDATDNS 

A  method  is  demonstrated  for  calculating  the  steady  heat  transfer 
rate  from  finned  cylinders  dissipating  energy  to  the  surroundings  by  con¬ 
vection  and  radiation.  Application  of  the  method  to  the  XM144E1  mortar 
revealed  that  the  radia+<ve  mode  accounts  for  60-701  of  the  total  energy 
transfer  under  some  operating  conditions.  Lack  of  adequate  convective 
heat  transfer  data  under  these  conditions  renders  the  calculations  somewhat 
tentative,  however,  and  it  is  recomnended  that  more  precise  experimental 
data  be  obtained  in  this  area. 

It  was  also  found  that  changing  the  fin  configuration  on  the  XM144E1 
mortar,  while  maintaining  the  same  weight  per  unit  tube  length,  can  further 
augment  the  heat  transfer  rate.  Strength  considerations  can  limit  the  gains 
achieved  through  changes  in  geometry,  however,  so  it  is  further  recommended 
that  appropriate  experimental  data  be  acquired,  probably  by  photoelastic 
methods,  for  predicting  the  load  carrying  characteristics  of  arbitrary  fin 
configurations.  This  data  can  then  be  used  as  a  guide  in  the  selection  of 
the  optimum  fin  configuration  —  optimum  from  the  viewpoint  of  heat  transfer. 
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APPENDIX  A  —  ANCLE  FACDRS 


The  following  angle  factors  were  obtained  using  the  standard  techniques 
given  in  ref.  [6]. 


*bl(Tb»  CO  "  1  (flitan<K  - 
w  /  2NR 


NR(aJ-b{) 
4ClCalcl+alel*dlbl)  "  a]tan«< 


(dibi-ajCiJai-eibi 


sin 


al  "  bl 


Nr 


(aj^bjcosej  ) 

*3. 

M  -  bf 


tan 
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tan^-  ^ 


-1  r 
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T1  "  ^  ’  (1’ 

f  1  *  NR  *  (trTb)2 
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•l  '  e>  Nr“"<  •  NR(’l‘Tb> 
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APPENDIX  B  —  NUMERICAL  METHODS 


General .  The  technique  for  determining  the  distributions  of 
(^(t^  end  Tj(^)  is  to  approximate  the  integrals  in  eqs.  (16) ,  (17)  end 
(19)  by  numerical  quadrature  formulas  and  then  solve  the  resulting  set  of 
non-linear  algebraic  equations  by  iteration.  The  method  employed  here  is 
basically  that  given  by  Donovan  (8]  and  the  reader  is  referred  to  his  work 
for  a  more  complete  discussion  of  the  technique. 


Sisson* s  one-third  rule  is  used  wherever  possible  in  the  form: 


r*in2 

F(>*)d*.  a  A*  2_  F(^i)c(i) 

i  ■  n. 


where  nj  and  n2  are  odd  numbers  with  n2>nj  ♦  2;  that  is,  Simpson's  rule  is 
applicable  only  to  a  set  of  double  panels,  each  panel  having  a  width  A .96. 
Also, 

"  b  ’  •  (■*) 
n2-  nx 

^  -  a  ♦  A^  (i-Oj)  nj<  isn2  (B3) 


i  -  n1# 
i  even 


(B4) 


In  equation  (Bl)  the  continuous  function  F(^)  is  replaced  by  a  set  of 
discrete  points,  F(^^),  uniformly  spaced  over  the  integration  interval  (b-a). 
In  the  present  application  the  fin  and  base  surfaces  are  subdivided  into 
(N-l)  and  (M-l)  intervals,  respectively,  as  shown  in  Figure  8.  The  width 
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of  each  interval  on  the  fin  is 


-  1  -  Nr 

N  -  1 

and  on  the  base 

dTb  »  [Nl  -  (l-Ng)tan  *(  ]  -  (1  -  NR)tan«< 

M  -  1 

■  NL  -  2(1-Nr)  tan  * 

M  -  1 

Both  N  and  M  Bust  be  odd  nuabers. 


CBS) 


(B6) 


Nuswrical  Approxiaation  of  the  Radiation  Equations,  Consider  the  application 
of  eq.  (Bl)  to  the  integral  in  eq.  (16a).  It  becomes 
1 

£  U[T*4(0  -T*4]  -  (l-eJQjttWjjllV  0<U 
R  J  ■  N 

(tJ4(J)-T*4]  -(l-e)Q1(J))K12(l.J)c(J)  (B7) 

J  -  1 

*  • 

where  T^(J)  and  Qj(J)  are  the  values  of  Tj(^j)  and  Q](fj)  *t  the  points 


^(J)  •  Nr  ♦  (J-l)  J  -  1,  2,  3,  ...»  N 

The  coefficients  c(J)  are  given  by  (B4)  with  nj  ■  1  and  n2  ■  N. 

The  function  K12(1,J)  refers  to  the  exchange  between  the  point  I  ■  1 
on  fin  surface  "1"  in  Figure  8  and  the  set  of  points  J  ■  lf  2,  3,  . N  on 
fin  surface  "2".  This  function  is  plotted  in  Figure  9  with  N  ■  21.  The 
curves  for  I  ■  3  and  5  are  also  shown.  The  values  of  Nj^,  and  •(  correspond 
to  the  configuration  shown  in  Figure  1. 
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It  is  noted  thst  Kj2(I»  J)  varies  quits  rapidly  for  snail  values  of 
I  and  J.  In  fact,  for  I  ■  1  the  peak  in  the  curve  is  not  well  defined  by 
the  darkened  circles  corresponding  to  the  points  J  ■  1,  2  and  3.  This  can 
result  in  significant  error  in  the  approxination  given  in  eq.  (B7).  The 
accuracy  of  the  approximation  can  be  improved  by  subdividing  the  two-panel 
integration  interval  1*J*3  into  the  eight-panel  interval  lsjs9  shown  in 
Figure  8.  The  angle  factor  is  designated  as  j)  and  is  plotted  with 

the  open  circles  in  Figure  9  for  I  ■  1  and  3.  This  function  is,  of  course, 
also  calculated  from  eq.  (A2). 

The  remainder  of  the  integrand  is  a  more  slowly  varying  function  of 
and  can  be  satisfactorily  approximated  at  the  points  j  ■  1,  2,  . ..,  9 
using  the  following  three-point  Lagrange  interpolation  formula: 

S(j)  •  p(p-l)  f(l)  ♦  (l-p2)f(2)  ♦  p(p*l)f(3)  (B8) 

2  2 


where 


P  -  (j  -  S)/4 


i  »  1*  2,  ...,  9 


and 


f(J)  -  c[Tj4(J)  -  T*4]  -  (l-cJQ^J)  J  •  1,2,3 


Using  this  notation  the  integral  in  eq.  (B7)  becomes 
1 

U[Ti4(C)  -  T*4]  -  Cl-e)Q1CO)K12(NR,OdC 


c 


«  A 


j  -  9 

S(j)KK12(l,  j)c(j)  ♦ 
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(B9) 


J  -  N 

:F«-! 


M(J)  -T*4]-(1-c)Q1(J))K12(1,J)c(J) 


The  first  summation  in  eq.  (B9)  extends  the  integration  over  eight  panels 
of  width  (A^/4),  that  is,  fro*  £  ■  NR  up  to  £  ■  Nr+  2A^.  The  second 
summation  covers  (N  -  3)  panels  of  width  Af  and  completes  the  integration 
fro*  C  *  ♦  2A^  up  to  £  ■  1.  Note  that  in  this  case  nj  »  3  and  n2  ■  N. 

All  of  the  integrals  in  eqs.  (16)  and  (17)  are  treated  in  the  sane 
Banner,  including  the  special  handling  of  the  integration  near  the  "corners 
The  complete  set  of  equations  used  to  determine  the  functions  Qi(^*i)  and 
Q^(t^)  are  as  follows: 


When  I  ■  1: 


j  -  9 


Q,tt>  -  e[l-Te4]  -  A£  ^S(j)KK12(l.j)c(j) 

4  j  ■  1 
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When  I  ■  4,  . ...  N: 


Qj (I)  -  c[T*4(I)  -T*4] 


J  -  N 


-  A(^T£[T*4(J)  -T*4]  -  (I-cJQj (J)  }K12 (I.J)c(J) 
J  -  1 


K  -  M 
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Ci  CD  l— 


K  ■  1 


When  K  •  1: 

%0)  -  e[l-T*4]  -  YsCDKKhjCUJcd) 

i  -  1 


I  -  N 


-  U[T|4(I)  -T*4]  -  (l-c)Q1(I)}Kbl(l,I)c(I) 
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(Bll) 


(B12) 


(B13) 
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In  eq.  (Bll)  T(k)  is  given  by 

TOO  ■  P(P-1)  g  (M-2)  ♦  (l-p2)*(M-l)  ♦p(p*l)g(M) 

2  2 

where 

P  *  SLZ  5)  km  i,2 . 9 

4 

g(K)  -  c (1-T*4)  -  (l-eJQjjOC)  K  ■  M-2,  M-l,  M 

The  functions  HX^^Ck,!)  and  KK^CK,!)  are  calculated  fro»  eq.  (Al)  at  the 
positions  indicated  in  Figure  8. 
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Numerical  Approximation  of  the  Integrated  Conduction  Equation.  Tha  two 
indafinita  intagrals  in  tha  integrated  conduction  aquation,  aq.  (19) 
involva  intagrating  ovar  an  odd  number  of  panels  of  width  A^  whan  I  is 
an  tvan  nuaber  and  over  an  even  nuaber  of  panels  whan  I  is  an  odd  nuaber. 
Thus,  Siapson's  rule  can  be  applied  directly  to  the  points  Tj(I), 

I  •  1,3,5,. . .  ,N,  but  an  alternate  quadrature  foraula,  applicable  to  a 
single  panel  of  width  A^  ,  Bust  be  employed  at  tha  even  numbered  points. 

Donovan  [8]  uses  the  following  foraula  froa  Haiaing  [11]: 

F(f)d/ «(A7«./24)[9FUi)  ♦  -  5F(>i*2A^) 

"+i 

♦  F(/*±*3A<jfr)] 


(B17) 


Specializing  Siapson's  rule  to  a  single  interval  of  width  2A<^  gives 
♦  2A* 

F(*)d*  MA^/SMF^)  ♦  4F(*^A.,0  ♦  F(^^2A^)] 


(B18) 


TB(Cl)  ’  A 

7N, 


and  applying  eqs.  (B17)  and  (B18),  tha  conduction  equation  becomes 
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T*(I)  -  l  ^[NR- fjCUlVi)  ♦  TB(I)  ♦  nr  ♦  [N^CjCDia^c  ) 
PCI)  (  6*  +  (l-NR)tan-(J 


(B19) 


where 


Ta(N)  -  0 

Ta(N-1)  -  (A(/24)[9*(N)  ♦  19*(N-1)  -S*(N-2)  ^(N-3)J 

When  I  »  N-2 ,  , 5,  3,  1 

Ta(I)  »  Ta(I+2)  ♦  (A^/3)  [♦(!)  ♦  4i|i  (!♦!)  ♦  4(1+2)] 


When  I  -  2,  4,  6,  ....  N-3 

ta(D  -  Ta(I-1)  -  (4^/24)  [9*(I-1)  ♦  194,(1)  -  5*(I+1)  ♦  *(I*2)] 

Also, 

tb(1)  -  0 

Tb(2)  -  (4^/24)  (94(1)  ♦  194(2)  -  54(3)  -  4(4)] 

When  I  ■  3,  5,  7,  . .«,  N 

TB(I)  -  TB(I-2)  ♦  (4^/3)  (4(1)  ♦  44(1-1)  ♦  4(1-2)] 

When  I  ■  N-l,  . ...  8,  6,  4 

V1)  *  VI+1)  *  (A{/24)[94(Iel)  ♦  194(1)  -  54(1-1)  +4(1-2) ] 

Note  that  Tj(l)  ■  1  and  need  not  be  evaluated  from  eq„  (B19). 

Method  of  Solution.  Eqs.  (BIO)  -  (B16)  and  eqs.  (B19)  represent  a  set  of 
simultaneous  non-linear  algebraic  equations  which  are  solved  by  the  method 
of  iteration.  A  weighted  average  of  the  two  most  recent  iterates  was  used 
to  accelerate  the  convergence.  This  took  the  form 
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Yi  ♦  2  •  °*6  4  °*4  n  ♦  i 

where  yj  4  j  and  y^  are  the  two  most  recent  iterated  values  of  an  tatknown 
and  Yi  +  2  represents  the  new  estiaate. 

For  the  fin  configurations  studied  in  this  report  convergence  was 
usually  achieved  in  less  than  fifteen  iterations.  It  should  be  noted, 
however,  that  for  very  tall  fins  (saall  values  of  NR)  or  for  large  tempera¬ 
ture  drops  along  the  fins  (large  N£  and  h#)  the  systea  aay  converge  very 
slowly  or  aay  not  converge  at  all.  This  difficulty  can  be  circumvented 
by  reducing  the  values  of  N£  and  h*  until  a  solution  is  obtained.  This 
solution,  along  with  slightly  larger  values  of  N£  and  h*,  is  then  fed  back 
into  the  program  as  a  "first  guess"  at  the  new  solution.  If  convergence 
is  achieved  the  process  is  repeated  until  the  desired  levels  of  N£  and  h* 
are  reached. 

In  comparing  the  computed  solutions  with  known  solutions  of  convective 
and  radiative  heat  transfer  from  fin  systeas  [1,2]  best  agreement  was  found 
when  the  step  sizes  Af  and  At^  were  set  approximately  equal.  Using  eqs. 
(B5)  and  (B6)  a  guide  for  choosing  N  and  M  is  then 
M  -  1  ^  Nl  -  2(l-NR)tanK 

N  -  1  1  -  nr 

Usually,  a  minimum  of  eleven  points  was  used  for  either  M  or  N  although 
smaller  values  can  and  should  be  used  if  the  above  guide  so  indicates. 

The  integrals  required  in  the  evaluation  of  the  fin  system  effective¬ 
ness  in  eq.  (24)  were  found  by  direct  application  of  Simpson's  rule,  eq. 
(Bl) .  The  calculation  of  the  actual  heat  transfer  rate  from  eq.  (25)  is 
straightforward.  A  listing  of  the  computer  program  is  given  in  Appendix  C. 
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APPENDIX  C  —  COMPUTER  PROGRAM 


A  litt  of  tho  Major  program  symbols  is  given  below.  The  symbol  ACC 

is  the  accuracy  criterion  employed  in  the  iteration  process.  Each  time 

* 

a  new  value  of  one  of  the  unknowns  Qj (I) ,  Tj(I)  or  (^(K)  is  calculated 
its  magnitude  is  compared  with  its  previous  value.  If  the  absolute 
value  of  the  difference  is  greater  than  the  accuracy  criterion,  ACC,  for 
any  one  of  the  unknowns,  a  further  iteration  is  calculated.  In  this  study 
the  value  ACC  »  0.001  was  used.  The  fin  angle,  ALPHA,  must  be  in  radians. 
Symbol  Used  in  Program  Symbol  Used  in  Text 


ACC 

ALPHA 

ATB 

C 

COND 

DELRO 

DELS 

DELTB 

EPS 

ETA 

HKB1(K,I) 

HS 

K12(I,J) 
KB1 (K, I) 


ACC 

< 

Tb 

C 

6# 

ATb 

c 

i> 

H*S>i(k»I) 

h* 

X12(I,J) 

SiCM) 


Sybol  Us«d  in  ProgT— 
KK12(I,J) 
KKB1(K,I) 

M 
N 
NC 
NL 
NR 
P(I) 

PHI (I) 

QKI) 

QB(K) 

QC 
QO 
QR 
QT 
RO 

ROl(I) 

S(J) 

T(K) 

TES 
T1S(I) 

X(I) 


Sybol  Us»d  in  T«xt 

kk12(ij) 

KKbl(,C'i) 

M 

N 


NR 

P(I) 

♦(I) 

Ql(l) 

Qc 

Qo 

Qr 

qT 

Ro 

fid) 

s(j) 

T(k) 

* 

Te 

TJ(I) 

*(I) 
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FLOW  CHART  OF  COMPUTER  PROGRAM 


Begin 


Read  in  and  print  out  input  data 


Set  initial  diatributlons  of  temperature  and  radiative  fluxes 


Calculate  new  fin  temperature  distribution.  Compare  each  new 
temperature  with  its  previous  value  to  see  if  it  satisfies 
accuracy  criterion;  if  not,  call  for  another  iteration  by  setting 
IDD  ■  1.  Also,  check  to  see  if  Tj(I)  is  greater  than  unity;  if 
sg,  set  it  equal  to  unity.  If  it  is  less  than  T*  set  it  equal  to 
Te.  (These  last  two  events  can  occur  in  an  iterative  scheme, 
especially  if  the  initial  guess  is  wide  of  the  mark.  By  restrict' 
ing  the  T|(i)  by  these  physical  bounds,  convergence  is  often 
achieved  in  a  system  that  might  otherwise  diverge.) 


alculate  base  radiative  flux  distribution.  Compare  each 
new  value  with  the  previous  value  and  call  for  another 
iteration  if  necessary. 


ux  distribution  and  check  for 
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